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Structured, periodic optical materials can be used to form photonic crystals capable of dispersing, routing, and 
trapping light. A similar phenomena in periodic elastic structures can be used to manipulate mechanical vi- 
brations. Here we present the design and experimental realization of strongly coupled optical and mechanical 
modes in a planar, periodic nanostructure on a silicon chip. 200-Terahertz photons are co-localized with me- 
chanical modes of Gigahertz frequency and 100-femtogram mass. The effective coupling length, Zom> which 
describes the strength of the photon-phonon interaction, is as small as 2.9 pm, which, together with minute oscil- 
lator mass, allows all-optical actuation and transduction of nanomechanical motion with near quantum-limited 
sensitivity. Optomechanical crystals have many potential applications, from RF-over-optical communication to 
the study of quantum effects in mesoscopic mechanical systems. 

Periodicity in materials yields interesting and useful phenomena. Applied to the propagation of light, periodicity gives rise to 
photonic crystals| 1 1, which can be precisely engineered to, among other things, transport and control the dispersion of light|2][3), 
tightly confine and trap light resonantly [4], and enhance nonlinear optical interactions[5|. Photonic crystals can also be formed 
into planar lightwave circuits for the integration of optical and electrical microsystems[6|. Periodicity applied to mechani- 
cal vibrations yields phononic crystals, which harness mechanical vibrations in a similar manner to optical waves in photonic 
crystals 1 7 8 9 [10] [TT] [12] [T3J . As has been demonstrated in studies of Raman scattering in epitaxially grown vertical cavity 
structures lfPfl and photonic crystal fibers] 15], the simultaneous confinement of mechanical and optical modes in periodic struc- 
tures can lead to greatly enhanced light-matter interactions. A logical next step is thus to create planar circuits that act as both 
photonic and phononic crystals| 16]: optomechanical crystals. In this spirit, we describe the design, fabrication, and charac- 
terization of a planar, silicon-chip-based optomechanical crystal capable of co-localizing and strongly coupling 200 Terahertz 
photons and 2 Gigahertz phonons. These planar optomechanical crystals bring the powerful techniques of optics and photonic 
crystals to bear on phononic crystals, providing exquisitely sensitive (near quantum-limited), optical measurements of mechan- 
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FIG. 1: Optomechanical crystal design, a, Geometry of nanobeam structure, b, Optical and c, mechanical bands and defect modes calculated 
via FEM for the projection of the experimentally-fabricated silicon nanobeam (A = 362 nm, w = 1396 nm, h y = 992 nm, h x = 190 nm, and 
; = 220 nm; isotropic Young's modulus of 168.5 GPa; n =3.493). In this particular structure, which will be referred to as "Device 1", 
^defect = 15 holes, ./V tota i = 75 and the spacing between the holes varies quadratically from the lattice constant of the projection (362 nm) to 
85% of that value (a "15% defect") for the two holes straddling the central cross-bar (the other parameters of Device 1 are as listed above). 
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ical vibrations, while simultaneously providing strong non-linear interactions for optics in a large and technologically-relevant 
range of frequencies. 

The geometry of the optomechanical crystal structure considered here is shown in Fig. [ij. The effectively one-dimensional 
(ID) optomechanical crystal consists of a silicon nanobeam (thickness t not shown) with rectangular holes and thin cross-bars 
connected on both sides by thin rails (we will refer to infinitely periodic constructs such as this as the "projection" of the 
finite structure). Fig. [lj> shows a finite-element-method (FEM) simulation of the optical band structure of the projection of a 
nanobeam (see caption for parameters). The electric field profile for modes at the band edge (k x = Jt/A, the boundary of first 
Brillouin zone) are shown to the right of the band structure. The finite structure terminates at its supports on both ends, forming 
a doubly-clamped beam. To form localized resonances in the center of the structure, the discrete translational symmetry of the 
patterned beam is intentionally disrupted by a "defect", consisting of a quadratic decrease in the lattice constant, A, symmetric 
about the center of the beam for some odd number of holes, Affect < A'totai- The defect forms an effective potential for optical 
modes at the band edges, with the spatial dependence of the effective potential closely following the spatial properties of the 
defect[ 17 1 (as illustrated in the inset of the optical band diagram). Thus the optical modes of the infinitely-periodic structure are 
confined by a quasi-harmonic potential. This effective potential localizes a "ladder" of modes with Hermite-Gauss envelopes, 
analogous to the modes of the ID harmonic potential of quantum mechanics. The localized optical modes of the finite structure 
(hereafter referred to as Device 1) are also found by FEM simulation and shown in Fig. [TJi to the right of the corresponding 
mode of the projection. 

Analogously, Fig. [TJ; shows a FEM simulation of the mechanical band structure of the nanobeam's projection. Mechanical 
modes at the band edge experience an effective potential analogous to the optical modes, localizing certain types of vibrations to 
the defect region. The colored bands give rise to mechanical modes that, when localized by the defect, yield "ladders" of modes 
with strong dispersive coupling to the localized optical modes (the frequency of the fundamental defect mode is indicated by 
a horizontal bar of the same color). The red mechanical band is analogous to the acoustic (longitudinal) vibrations of a solid 
or system of masses and springs, starting off at k x — as a pure translation of the entire structure and evolving to differential 
acoustic vibrations of every pair of cross-beams at the band edge (mode at k x =7t/A shown to right of band diagram), just as 
the highest-frequency acoustic vibrations of a solid have every unit cell vibrating out of phase with its nearest neighbors. The 
green band also involves in-plane vibrations of the cross beams, but, in this case, the rails recoil in opposition to the cross- 
beam vibration, which means that the structure can vibrate (rather than translate) even at very long-wavelengths, giving rise to 
a non-zero frequency at k x = 0. The blue band consists of in-plane vibrations transverse to the axis of the structure, where, at 
r (mode at Y shown to right) the structure appears to "breathe", as the cross-beams are stretched and compressed. We classify 
these optomechanically-coupled mechanical modes, from lowest to highest frequency, as "pinch", "accordian", and "breathing" 
modes. The localized mechanical modes of Device 1 are shown to the right of the corresponding mode of the projection. 

The two kinds of waves, mechanical and optical, are on equal footing in this structure. Each mechanical mode has a frequency 
V m = £l m /2% and displacement profile Q(r); each optical mode has a frequency v = C0 o /27i and electric field profile E(r). 

Just as the optical mode volume, V = / dV ^ — ^/gE|) ^ ' describes the electromagnetic localization of the optical mode, 

the mechanical mode volume, V m = p J dV ( ma j^Q|) ) ( see App. Di, describes the strain energy-averaged localization of the 
mechanical mode. For both the localized optical and mechanical modes of the patterned beam cavity, the effective mode volume 
is less than a cubic wavelength. The effective motional mass, being proportional to the mode volume (m e ff = pV m ), is between 
50 and 1000 femtograms for the mechanical modes shown in Fig. [Tp (ps, = 2.33 g/cm 3 ). 

Drawing on recent work in the field of cavity optomechanics lfT8lll9L we describe the coupling between optical and mechanical 
degrees of freedom (to lowest order) by an effective coupling length Lom = ( ^ ) 1 ( see App. pL where 8v is the change in 
the frequency of an optical resonance caused by the mechanical displacement parameterized by aTTor this work, a is defined as 
the maximum displacement that occurs anywhere for the mechanical mode. By definition then, the smaller Lom, the larger the 
optical response for a given mechanical displacement. Lom is also the length over which a photon's momentum is transferred 
into the mechanical mode as it propagates within the structure, and thus is inversely proportional to the force per-photon applied 
to the mechanical system. 

To calculate Lom, we employ a perturbative theory of Maxwells equations with respect to shifting material boundaries 1 20 1. 
The derivative ^ around some nominal position, where the optical fields are known, can be calculated exactly without actually 
deforming the structure for a surface-normal displacement of the boundaries, h(a;r) = Q(r) n = ocq(r) • n, where q(r) = 
Q(r)/a = dQ(r)/da is the unitless displacement profile of the mechanical mode, and a parameterizes the amplitude of the 
displacement. Using this perturbative formulation of Maxwell's equations, we find 
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FIG. 2: Photonic and phoninic crystal mode spectroscopy, a, and b, show SEM images of the fabriced silicon nanobeam optomechanical 
crystal, c, Optical spectroscopy of Device 1 with the taper waveguide in contact, d, Mechanical spectroscopy of Device 1 with taper waveguide 
in contact, e-g, Zoomed RF mechanical spectra of Device 1 showing pinch (red), accordian (green), and breathing (blue) modes, h, Funda- 
mental breathing mode RF spectra (of a third device) at low (15 /jW) and high (190 /jw) coupled optical power . At high optical power the 
breathing mode is optically amplified through dynamical back-action of the radiation pressure force, with an above threshold effective g-factor 
greater than 10 . 



where n is the unit normal vector on the surface of the unperturbed cavity, D(r) = e(r)E(r), Ae = £j — £?, A(e~') = £j — , 
£i is the dielectric constant of the periodic structure, and £2 is the dielectric constants of the surrounding medium (£2 = £0 in this 
case). This method of calculating the coupling provides a wealth of intuition about the nature of the coupling and can be used to 
engineer the structure for strong optomechanical coupling. 

Figs. [2^ and |2j? show scanning electron microscope (SEM) images of a fabricated silicon nanobeam with the parameters of 
Device 1 (see App. [Aj. The optical modes of the nanobeam are probed with a tapered and dimpled optical fiber[21 1 in the near- 
field of the defect cavity, simultaneously sourcing the cavity field and collecting the transmitted light in a single channel (see 
App. [B]l. Fig. [3J: shows the low-pass filtered optical transmission spectrum of Device 1 at low optical input power (~ 30 /jW). 
The optical cavity resonances are identified by comparison to FEM modeling of the optical modes of the structure. Looking 
in the radio frequency (RF) spectrum provides information about the mechanical modes of the structure, as mechanical motion 
gives rise to phase and amplitude modulation of the transmitted light. Figs. [2j;-f show the measured photodetector RF power 
spectral density (PSD) of the optical transmission through the second order cavity resonance (this mode was used due to its 
deep on-resonance coupling). A series of lower frequency modes can be seen in the spectra (~ 200 MHz and harmonics), 
corresponding to compression of modes of the entire beam, followed by groups of localized phononic modes of the lattice at 
850 MHz (pinch), 1.41 GHz (accordian), and 2.25 GHz (breathing). The transduced signal at low optical power corresponds to 
thermally-excited motion of the mechanical modes, and is inversely proportional to m e ffL,Q M (see below). At higher optical input 
power (see Fig.[2ji), optical excitation of regenerative mechanical oscillation! 18 1 of the breathing modes is possible due to the 
small mass and short optomechanical coupling length of the co-localized phonon and photon modes. 

Fig. [3^, top panel, shows how the frequency of the fundamental mechanical breathing mode scales with a uniform geometric 
scaling in the plane. A series of 12 devices have been fabricated, identical except that the entire geometry in the plane is 
scaled incrementally by 1% per device. For each device, one of the first two optical modes is selected and used to measure the 
mechanical frequency of the fundamental breathing mode (cyan dots, top panel). The frequencies plotted in magenta are the 
normalized frequencies, i.e. the bare frequencies (cyan) times the scale factor for the device. Because the vibrations are entirely 
two-dimensional (in the plane), the frequency of the mechanical mode scales perfectly with the two-dimensional scale factor. 
This is in contrast to the optical modes (Fig. [3^, bottom panel), which clearly do not scale with the planar geometry, a result 
of the coupling of in-plane and vertical optical mode confinement (scaling in all three dimensions is thus required). Since the 
planar scaling for the lattice-localized mechanical modes is trivial, this method could be used with a larger span of devices to 
measure the frequency dependence of the Young's (or bulk) modulus of the material. 

Significant shifts in the frequency of the lattice-localized mechanical modes can be obtained through a non-uniform planar 
scaling. Fig. [3jB shows the RF PSD for Device 1 and a second device, Decice 2, which has an essentially identical lattice 
constant, A = 365 nm, and total length, L, as compared to Device 1, but a considerably smaller width (w — 864 nm, h y = 575 
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FIG. 3: Phononic mode tuning and transduction, a, Geometric scaling (planar) of the fundamental breathing mode. Device 1 is the 
device with scale factor 1.03. The best linear least-squares fit lines in the top panel correspond to the mechanical frequency changing by 
—0.9% ± 0.2% per device; the normalized frequency changes by —0.01 % ± 0.2% per device. The optical frequency of the mode used to make 
the mechanical measurement is filled (the other optical mode is open), b, Engineering of pinch mode frequencies, showing two devices with 
pinch mode frequencies of 850 MHz and 1.75 GHz. The mechanical band diagrams of each structure are shown to the right of the measured 
RF spectrum, with the pinch mode band highlighted in red. c, Transduction of breathing mode motion, d, Ideal and actual ("primed") modes 
of the silicon nanobeam optomechanical crystal due to the ideal (dashed) and actual defect (solid, red). 



nm, h x =183 nm). Simulations show that the pinch modes are the lowest-frequency group of localized and optomechanically- 
coupled mechanical modes in both structures (see right panels of Fig. |3j3). Experimentally, the ratio of the localized pinch mode 
frequencies (highlighted in red) in these two devices is 1.749 GHz/805 MHz = 2.17. The ratio of the frequency of the localized 
pinch-mode manifold, after accounting for the defect, is theoretically 1.826 GHz/846 MHz = 2.16. It is interesting to note that 
the mechanical modes of the entire doubly-clamped beam (as opposed to the lattice-localized modes) depend very weakly on the 
structural differences betweenJDevice 1 and Device 2. For instance, the second-order acoustic vibration mode of the nanobeam 

(highlighted in yellow in Fig. 



3 3) has a frequnecy which should be ^ J ^y, where E is Young's modulus and (p) is the average 

linear density. The frequency of this mode is measured to be 234 MHz/195 MHz = 1.20 times higher in Device 2 than for Device 
1, which is in good agreement with the ratio -\/(pi) / (P2) = 1.23. The difference between the change in the frequencies of 
the lattice-localized versus beam modes illustratres the independence of these two "systems"; once the wavelength of the global 
beam modes approach the scale of the lattice periodicity, the vibrations become localized and behave independently of the global 
beam structure (such as the end clamps). 

Fig. [3}; shows the RF optical transmission spectrum due to Brownian motion of the breathing modes of Device 1 (i.e., at low 
optical input power), for the three optical modes shown in Figs. [TJ; and Because the various optical modes have different 
spatial profiles, each mechanical mode has a different Lom far each optical mode. The root-mean-square (rms) mechanical 
amplitude of a mode due to Brownian motion is (a 2 ) = ksT / (m e ffQ. 2 ) . It can be shown analytically that the factor 1 / (m e ffLo M ) 
uniquely determines the transduction of the Brownian motion for these sideband-resolved optomechanical oscillations (see App. 
Je]i. To the right of each measured spectrum is the experimentally-extracted mechanical frequency and value of l/(m e ffLQ M ), 



together with the values of these quantities obtained from the FEM model (using Eq. D4 to determine Lom)- Good corre 



spondence, in both frequency and transduced signal amplitude, is found across all optical and mechanical mode pairs. In order 
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to achieve this level of correspondence, imperfections in the fabricated structure are taken into account by extracting the ge- 
ometry from high-resolution SEM images of the device and calculating the modified optical and mechanical modes (Fig. [3}i). 
The resulting measured value for optomechanical coupling between the fundamental breathing and optical mode (assumming 
a FEM-calculated motional mass of m e ff = 330 fg) is Lom = 2.9 //m, approaching the limit of the wavelength of light. The 
sensitivity of the mechanical transduction of the fundamental breathing mode can be appreciated by comparing the mode's rms 
thermal amplitude at T = 300 K, a,/, = 245 fm, to its quantum zero-point motion of a zp = 3.2 fm. The sensitivity limit, as given 
by the background level in the middle panel of Fig.|3j;, is thus a factor of ~ 7.5 times that of the standard quantum limit. 

The loss of mechanical energy from confined mechanical modes of a phononic crystal can, in principal, be made arbitrarily 
low (and thus the mechanical Q arbitrarily high) by including a large number of unit cells outside the localizing potential region. 
Of course, other forms of mechanical damping, such as thermo-elastic damping, phonon-phonon scattering, or surface damping 
effects, would eventually become dominant|22|. This makes optomechanical crystals ideal structures for studying these loss 
mechanisms. The fundamental breathing mode of the ID phononic crystal structure studied here, at 2.254 GHz, has a room 
temperature mechanical Q of 1300 in air, and in contact with the taper waveguide (power-dependent measurements confirm that 
this mechanical Q is not enhanced by dynamical back-action). This corresponds to a frequency-^* product of 3 x 10 12 Hz, a value 
close to largest demonstrated to date [ 23 1. Although further tests (as a function of temperature and lattice periods) are required 
to determine the contribution of various mechanical loss mechanisms, numerical simulations show that mode coupling between 
localized and leaky phonon modes exist in these ID cavity structures and can significantly limit the g-factor (see App. |Cj. This 
obstacle can be overcome in two-dimensional periodic slab structures, which have been shown to possess complete gaps for both 
optical and mechanical modes simultaneously J24]. 

The experimental demonstration of optomechanical coupling between 200 Terahertz photons and 2 Gigahertz phonons in a 
planar optomechanical crystal paves the way for new methods of probing, manipulating, and stimulating linear and non-linear 
mechanical and optical interactions in a chip-scale platform. As the study of quantum mesoscale mechanical oscillators has 
nearly become a reality (251 l26l l27l l28l [29], high frequency mechanics will provide a distinct experimental advantage due to 
the lower thermal phonon occupancy. In addition, optomechanical crystals with full phononic bandgaps provide a platform to 
decouple the direct decoherence (phonon leakage) of mechanical modes from their supports. This could allow the preparation 
of mechanical vibrations with ultra-long lifetimes, the study of the intrinsic mechanical material losses, and narrow-linewidth 
Gigahertz frequency sources. Optomechanical crystals could also be used as high-spatial resolution mass sensors; with m e g = 
62 fg and v m = 850 MHz, the mass of a single Hemoglobin A protein (~ 10~ 19 g) would change the frequency of the pinch 
mode by 700 Hz, allowing sensitivity paralleling NEMS zeptogram mass sensors[30|. 
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APPENDIX A: FABRICATION 



The optomechanical crystal nanobeam is formed in the 220 nm thick silicon device layer of a [100] Silicon-On-Insulator (SOI) 
wafer. The pattern is defined in electron beam resist by electron beam lithography. The resist pattern is transfered to the device 
layer by an inductively-coupled plasma reactive ion etch with a C4Fg/SFg gas chemistry. The nanobeam is then undercut and 
released from the silica BOX layer by wet undercutting with hydrofluoric acid. 



APPENDIX B: EXPERIMENTAL SETUP 
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FIG. 4: Experimental setup used to measure optical, mechanical, and optomechanical properties of silicon optomechanical crystal nanobeam. 
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The experimental setup used to measure the optical, mechanical, and optomechanical properties of the silicon optomechanical 
crystal nanobeam is shown in Fig. [4] The setup consists of a bank of fiber-coupled tunable infrared lasers spanning approximately 
200 nm, centered around 1520 nm. After a variable optical attenuator (VOA) and fiber polarization controller (FPC), light enters 
the tapered and dimpled optical fiber, the position of which can be controlled with nanometer-scale precision (although vibrations 
and static electric forces limit the minimum stable spacing between the fiber and device to about 50 nm). The transmission from 
the fiber is (optionally) passed through another VOA and finally reaches an avalanche photodiode (APD) with a transimpedance 
gain of 11,000 and a bandwidth (3 dB rolloff point) of 1.2 GHz. The APD has an internal bias tee, and the RF voltage is 
connected to the 50 Ohm input impedance of the oscilloscope. The oscilloscope can perform a Fourier transform (FT) to yield 
the RF power spectral density (RF PSD). The RF PSD is calibrated using a frequency generator that outputs a variable frequency 
sinusoid with known power. 

APPENDIX C: SIMULATION PARAMETERS 

Modeling of both the optical and mechanical modes is done via finite element method (FEM), using COMSOL Multiphysics. 
The following subsections provide the description of the method used to define the FEM model of the optomechanical crystal 
system. 

1. Extracting the geometry in the plane 

To model the optomechanical crystal system, the geometry of the as-fabricated structure must be measured. As the features 
are smaller than an optical wavelength, the measurements must be done by scanning electron microscope (SEM). Figure [5^ 
shows an "eagle's-eye" high-resolution SEM micrograph of a portion of device 1, with the defect centered in the image. 



a 




FIG. 5: (a) Scanning electron micrograph of fabricated silicon optomechanical crystal, (b) Approximated geometry shown as blue overlay on 
SEM micrograph from (a). 



Digital line-scans of the micrograph are used to detect the edges of the geometry. From the extracted edge positions, the 
geometry is approximated as a series of rectangular holes with two filleted ends inside of a rectangle (the beam), giving an 
approximate planar geometric representation of the structure shown as an overlay in Fig.[5]3. This geometric representation takes 
into account the size, position, and any curvature of each hole, giving an accurate approximation of the geometry. In the defect 
region, each hole is given by its measured value. Outisde the defect region, a series of holes is used to get the average hole shape, 
which is used in the model. 

The SEM has been calibrated, and the dimesnions as measured by the SEM are too large by 5%. Thus, the entire planar 
geometry is uniformly scaled down by 5%. Since the lattice constant, A, is a center-to-center distance between features, it is not 
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affected by erosion during processing, which makes it the most reliable measure of distance on the sample. After applying the 
SEM calibration factor, the average lattice constant outside of the defect as measured by the SEM agrees with the value written 
by the electron beam lithography tool to better than 1%. Since the SEM and lithography tool are independent, this is yet another 
confirmation that the geometry has been measured correctly (the fine spectral features of the simulation are the other way to 
check the geometry measurements, after comparing to measured mechanical and optical spectra). 

The SOI wafer thickness is specified as 220 nm by the manufacturer. We will assume that the planar geometry extends 
uniformly into the vertical direction for the entire 220 nm, since the "eagle's-eye" view used to measure the planar geometry 
does not capture any asymmetries in the vertical dimension. These vertical asymmetries are much more difficult to extract 
without sacrificing the device (by focused ion beam, cleaving, etc.). 



2. Young's modulus and index of refraction 

The nanobeam structures are fabricated such that the long axis (i) is parallel to the SOI wafer flat, which is oriented along 
[110] (±0.5°). We decompose the displacement field in FEM simulations along the crystal axes and find the majority of the 
strain energy is primarily stored in deformations along the family of equivalent directions specified by (1 10). Because the strain 
for the modes of interest are primarily along (110), an isotropic elasticity tensor derived from a single Young's modulus and 
Poisson's ratio is appropriate for the current level of detail. The index of refraction will also be treated as an isotropic scalar. 

As the Young's modulus, E, and index of refraction, n, determine the phase velocity of the waves (and thus the frequency), 
they can be "tuned" to make a single simulated frequency (mechanical for E and optical for n) come out exactly as measured. 
The free spectral range of the modes and the relative frequencies of different types of modes are determined by the details of 
the geometry, in conjunction with E and n; so although a single frequency can always be made to match experiment exactly 
by scaling E or n, the wider details of the spectrum are a more accurate reflection of whether the model is a good match to the 
experimental values. 

After accounting for the planar geometry and scale factor, the Young's modulus and index are tuned until the fundamental 
optical mode and the fundamental breathing mechanical mode each come out exactly as measured, which occurs for E = 168.5 
GPa and n = 3.493. These parameters, along with the measured geometry (as discussed above) yields a model that produces the 
values shown in Fig. [3J;. 

3. Optics: mode maps and modeling 
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FIG. 6: (a) Optical modes measured in a 200 nm laser wavelength span for a series of 20 devices, (b) Simulation of Device 1, which is the 
device with scaling factor 1.03 (dashed line). Filled blue circles correspond to modes of the fundamental (valence) optical band; the region 
shaded blue corresponds to frequencies that are no longer within the defect potential barrier height (i.e. propagating modes). Open blue circles 
correspond to transverse valence band modes. Black circles correspond to conduction band modes. 
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Fig. [6^ shows all the optical modes measured for a series of 20 devices, which are identical up to a uniform planar scaling 
that changes by 1% per device (Device 1 is the device with scale factor 1.03). Because of the limited laser range, only a limited 
number of modes can be measured on any given device. By measuring this series of uniformly scaled devices, a large number 
of modes can be seen. 

The devices, taken together, display a number of conspicuous features. First, the devices display a series of five modes with 
relatively high optical Q (relative to the low-g, waveguide-like modes that appear outside this range). Second, the smallest 
devices show a number of low-g modes at frequencies below the fifth mode. Finally, at high frequencies, the devices display 
another set of low-g modes, which are higher-g than waveguide-like modes but not as high-g as the other five modes. 

These features are all consistent with the optical model of Device 1 . Figure[6j5 shows the simulated modes of Device 1 , plotted 
as a function of their optical Q. The simulation shows that the defect confines 5 modes, with a precipitous drop in Q as the modal 
frequencies exit the defect potential (go above the energy barrier height). The simulation also explains the series of modes higher 
in frequency, which are not the conduction band modes. These modes, indicated with open circles in both Fig. [6^ and [(J (as 
opposed to filled) circles, are the Hermite-Gauss ladder of modes with a single node transverse to the direction of propagation. 
These modes have a lower effective index, which reduces their radiation-limited Q relative to modes without transverse nodes. 
Conduction band modes (which are not measured) are shown as filled black circles. 

Even though the optical information provided by any single device would be difficult to unravel, the measurements of the 
series of devices coupled with simulation allow us to unambiguously identify the optical spectra of every device in the series. 

4. Modeling mechanical Q 




FIG. 7: (a) In-phase and (b) in-quadrature mechanical displacement field of fundamental breathing mode with absorbing "pad". 



Simulations of an isolated optomechanical crystal nanobeam (hard boundary conditions at the edges) show that the confined 
modes can couple to modes of longer wavelength, the effect of which is to produce a standing wave within the structure that 
is not attenuated outside the defect. In the real structures, these long-wavelength modes will travel down the nanobeam and 
partially reflect at the contacts due to an effective impedance mismatch between the nanobeam and the bulk, with the rest of 
the power leaving the structure into the bulk. Thus the isolated mode is coupled to a mode with long wavelength but identical 
frequency that can radiate into the surrounding "bath". 

To model the loss due to this resonant coupling to radiative modes, we include a large, semi-circular silicon "pad" on each side 
of the nanobeam. To make the pad act like a "bath", we introduce a phenomenological imaginary part of the speed of sound in 
pad region; i.e., v pat j — > vsiii CO n(l + /T|), where v = •v/E/p. This creates an imaginary part of the frequency, and the mechanical 
Q can be found by the relation, Q m = Re{v m } / (2Im{v m }). 

By adding loss to the pad material, part of the power reflected at the contacts will be due to the change in the impedance 
from the absorption. From this point of view, r| should be made as small as possible, since this contribution to the reflection 
coefficient is not present in the real system. However, r| must also be large enough that the self-consistent solution includes a 
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radiated wave that propagates for a significant portion of the pad, which only happens if the wave is appreciably attenuated by 
the time it reflects from the edge of the simulation and returns to the contact. Thus, the pad is made as large as possible, given 
computational constraints, and the absorption is increased until the Q m changes appreciably, which gives the threshold value for 
T| at which the reflectivity of the contacts has an appreciable contribution from the absorption. The simulation is thus performed 
with a value of r| that produces a propagating wave in the pad without causing an artificial reflectivity at the contact; propagation 
in the pad is easily verified if the position of the nodes/antinodes swap between the in-phase and in-quadrature parts of the 
mechanical cycle (the nodes/antinodes of a standing wave are stationary). Figs. [7^ and[7]3 show the in-phase and in-quadrature 
(respectively) parts of the optical cycle, showing a propagating radiative mode in the pad. The simulated values of Q m can be 
found in Table [I] Most of the values of mechanical Q calculated this way are in fair agreement with the measured values except 
for the fifth breathing mode, which, by this method, appears to have a Q m of over 22,000; in reality, the Q of this mode may be 
limited by coupling with a leaky mode caused by defects that are not modeled (such as vertical defects), material losses, or the 
presence of the taper waveguide. 

Table[I]summarizes the properties of the breathing mechanical modes. Measured values are denoted with a tilde and extracted 
using techniques described below. Simulated values are calculated using methods also described here in the appendices. 



TABLE I: Measured and Simulated properties of the breathing mechanical modes. Tildes indicate measured quantities. The experimental 
effective lengths, Lqm, between each breathing mode and the first three optical ca vity modes are calculated using the experimentally extracted 
m e ff^0M ( see FiS-pb of main text) and dividing by the »ieff from the model. See 5 C4 for discussion on modeling Q m . 



Mode 


v m (GHz) 


v^ (GHz) 


™eff (fg) 






2 ^OM (Mm) 


2 ^OM (<" m ) 


3 Lom 0"m) 


3 ^OM C"m) 


Qm 


Qm 


1 


2.254 


2.254 


329 


4.9 


2.9 


6.4 


5.1 


7.8 


4.4 


1280 


2050 


2 


2.275 


2.270 


399 


7.1 


4.5 


6.2 


12 


9.5 


9.8 


1130 


1180 


3 


2.294 


2.290 


628 


11 


N/A 


6.2 


5.3 


7.7 


4.1 


613 


1290 


4 


2.322 


2.326 


704 


110 


N/A 


64 


49 


26 


N/A 


973 


387 


5 


2.369 


2.361 


665 


38 


N/A 


11 


25 


7.1 


4.7 


950 


21600 



APPENDIX D: CAVITY OPTOMECHANICS: DEFINITIONS 



For clarity and self-consistency, the following sections will present all the aspects of cavity optomechanics relevant to the 
experiments, in the same notation used in the main text. Much of the information can be found in other places lfT~8l [3T1 [32l l33ll . 



1. Relationship Between Fields and Their Corresponding Profiles and Amplitudes 

Cavity optomechanics involves the mutual coupling of an optical mode to a mechanical mode of a deformable structure. 

The optical mode is characterized by a resonant frequency C0 o = 2nv and electric field E(r). The mechanical mode is 
characterized by a resonant frequency £l m = 2nv m and displacement field Q(r), where Q(r) is the vector displacement describing 
perpendicular displacements of the boundaries of volume elements. The cavity optomechanical interactions of the distributed 
structure and its spatially-dependent vector fields can be reduced to a description of two scalar mode amplitudes and their 
associated mode volumes, with the coupling of the amplitudes parameterized by a single coupling coefficient. 

The mode amplitude, c, and complex vector field profile, e(r), are defined such that the complex electric field is E(r) = ce(r) 
(the physical field is given by the real part of E(r)e ! ™)- For pedagogical reasons, the amplitude c is normalized such that the 
time averaged electromagnetic energy is equal to \c\ ; i.e. U = \c\ = \ J dVe |E| . This forces e to be normalized such that 1 = 

/ \ 2 

\ J dVe|e| 2 . In cavity quantum electrodynamics, one also defines an effective optical mode volume, V = J dV 



y£|E 



V 



in order to gauge the strength of light-matter interactions. 

The acoustic field amplitude, a, and mode profile, q(r), are defined such that Q(r) = ccq(r). a is defined as the largest 
displacement that occurs anywhere for the mechanical field, Q(r), is a. It is important to note that this particular choice of a 
determines the mechanical mode effective volume and effective mass, V m and m s s = pV m , respectively. In particular, this choice 

of a requires the complimentary definition: m e g = p / dV ( ma j^Q|) ) ■ To see this, note that the free evolution of the mechanical 

oscillator has, by definition, a time-independent total energy ^mechanical = "^(n 2 oc 2 + a 2 ). On the other hand, integrating the 
total energy of each volume element must also give this same total energy. If we pick the point in phase space at which all 
the mechanical energy is potential energy (i.e. the turn-around point), we must have that ^mechanical = J p|Q( r )| 2 dV = 
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jm e ff£2 2 a 2 , or, in other words, m e ff<X 2 = / p|Q(r)| 2 dV. One can arbitrarily choose the definition of the amplitude or the mass, 
but choosing one determines the other. For a system like a confined mode of a phononic crystal defect cavity, where only a very 
small portion of the total mass undergoes appreciable motion, the most sensible choice of the mass is the amplitude-squared 
weighted density integral, which, as stated above, is the choice of mass associated with a = max(|Q(r)|). 

2. Coupling Between the Mechanical and Optical Modes 

The optomechanical coupling affects the optical mode by tuning its resonant frequency as a function of displacement, C0 o (a); 
whereas the coupling affects the mechanical mode by applying a force, which is expressed as a gradient of the cavity energy, 
d \c\ 2 /da. The optical resonant frequency is usually expanded in orders of the (small) displacement, a around some equilibrium 
displacement, oco. 



«o(a) 



CO, 



+ (a-a ) — 
a=ao act 



a=ao 



+ (a-a ) 2 



d 2 (0 o 
da 2 



a=ao 



(Dl) 



In the case that the terms higher than first order can be neglected, this equation simplifies to 



coo(a) 



CO, 



+ (a-a ) — 
a=ao dot 



= co + (a - a )goM = co + (a - a ) - — 

a=a Lqm 



(D2) 



where co = co 



a=a 



is the equilibrium resonance frequency of the optical mode, goM 



da 



is the derivative of the 



a=a.Q 



a=ao 



is the effective length of the optome- 



resonance frequency of the optical mode evaluated at equilibrium, and Lq^ = i ^ 

chanical system, which is a universal parameter that relates the displacement relative to Lom to a relative change in optical 
frequency (i.e. oc/Lom = §0) /co ). It is simple to show that Lom is equal to the length of a Fabry-Perot cavity with one fixed 
mirror or the radius of a microtoroid or microdisk for the radial breathing mode. For a zipper cavity or double-microdisk, Lom 
is a function of the spacing between the coupled elements that is exponential when the two elements are in the near-field of each 
other, with Lom approaching the wavelength of light as the spacing approaches zero. 
The optical force is the gradient of the optical mode energy 



-optical | 



d|c| z _ d|c| z dco _ Ic^dcoo 
da dco da co da 



Lom 



(D3) 



The perturbation theory of Maxwell's equations with shifting material boundaries 1 20 1 allows one to calculate the derivative of 
the shift of the resonant frequency of the optical modes of a structure in terms of some parameterization of a surface deformation 
perpendicular to the surface of the structure. Thus, if the result of a mechanical simulation is the displacement field, Q(r) = 
aq(r)=aQ(r)/max(|Q|) 



1 

Lom 



(14 



dq 

da 



Ae|E|if-A(e" 1 )|D J 



(D4) 



dVelEl- 



where n is the unit normal vector on the surface of the unperturbed cavity, Ae = £i — £2, A(e 1 ) 
constant of the periodic structure, and £2 is the dielectric constants of the surrounding medium. 



£ 9 , £1 is the dielectric 



3. Coupled Equations of Cavity Optomechanics 



For an optical field E(r), which is described by the mode amplitude c, and a displacement (acoustic) field Q(r), described 
by mode amplitude a, the optomechanical interaction of linear order (terms of order d 2 co /da 2 and higher are neglected and 
ao = 0) of the fields is governed by the coupled differential equations 
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*(*)= (-f-f«)o(lT^))c(0+^ e - to (D5) 
a + r,a + ^ a=±w j!_ (D6) 

where K = K; +K ex t, K; = tt>o/Gop.i is the intrinsic optical loss rate of the cavity; Ke X t is the extrinsic coupling rate between 
the optical input field and the optical cavity field; s is the amplitude of the input field, normalized such that \s\ 2 = Pq is the 
optical input power impinging on the cavity; £2 is the acoustic resonance frequency of the cavity |[34l ; T = £l/Q m t _ is the intrinsic 
acoustic loss rate of the cavity; and m t _j = pV a coustic is the effective mass of the acoustic mode of the cavity, being the product of 
the mass density and the effective volume of the acoustic mode. Here we use the convention that the optical resonance frequency 
is decreasing with increasing a and the optical force is the positive gradient of the cavity energy; this is the typical convention, 
but it is completely arbitrary l35ll . 



APPENDIX E: OPTICAL TRANSDUCTION OF MECHANICAL MOTION 



One approximate method of solving the linear equations of optomechanics (D5 1 and (D6 1 is to treat the acoustic mode as a 
small perturbation of the optical mode, find the effect of the acoustic mode on the optical field, and then calculate what effect 
the perturbed optical mode has on the acoustic mode. Starting with the premise that the displacement, a, is sinusoidal; i.e. 



the mode amplitude, c, is thus described by 



a(t) = aosin(£2f) , (El) 



lcW+yf*-. (E2) 



1. Formal Solution 



The homogeneous solution to (|E2[) is 



(f)=Coexp^(-|-/co )f-/-^^cos(£2r)^) . (E3) 



The particular solution to (E2 1 can take the form c p (t) = C(f)c/,(f), satisfying 



1 • [*e -iw ■ i^e _____ ( , . \ , . 0C C0 o 



C(/) = — w se = is \ -tt exp - + icoo f + iy^- " cos(n?) - mt , (E4) 
c h {t) \ 2 I 2 \\2 ) Lqm £2 



where Co has been absorbed into C(f). 
Introducing the modulation index 

Lom £2 

the cosine part of the argument of the exponential can be expanded into Bessel functions as 

rap(±iPcos(flf))= £ (±i) n J n (V)e inat , (E6) 

n=—oo 

which is known as the Jacobi-Anger expansion. This allows straight-forward integration of C; to wit, 
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C(t) = is 



n=-~ 2 



fV»(P) 



+ i(nQ-A) 



,;'(f+«n-«+a>o)l 



(E7) 



Thus 



C p (t)=C(t)c h (t)=isJ^ £ lc~T~ 



i(nQ-G>)r+#cos(Qf) 



-i(nQ-A) 



(E8) 



The general solution is then c(f) = c/,(f) + c p (f). As c/,(f) is exponentially damped at rate K, the general solution rapidly 
converges to c p (t), which is the steady-state solution. This optical mode amplitude can thus be used to compute the various 
properties of the optomechanical system. 



2. The Transmission of an Oscillating Cavity 

The steady state power exiting the cavity is 



s out = ise ~~ - \/ —c p (t) 



(E9) 



and thus 



\ s out I — I ls out I 



se 



= N 2 +|M0| 2 -2Im 




(E10) 



with 



-21m 




c„s e 



(n—m) 




i(n+m)Q.t 

'A) 



(Ell) 



and 



/("- ffl )y„(p)y m (p) 



n,m 



(| + i(nQ. -A)) (| -/(mil -A)) 



J(n—m)Q.t 



(E12) 



3. RF Spectrum of the First Order Sidebands 



The experimental arrangement in this work is such that the optical mode is observed by weakly populating the cavity with 
photons via the tapered optical fiber waveguide and then collecting the transmitted photons via the same waveguide. The RF 
spectrum of the power transmitted, \s out \ 2 , is obtained by simple photomixing on an APD. This RF spectrum contains information 
about the mechanical modes because the mechanical modes modulate the optical mode via the optomechanical coupling. The 
modulation manifests as a set of sidebands created in the cavity that are transmitted and detected at the APD. If the amount of 
optical modulation, p, is small (i.e. p <C 1), only the first sideband contributes, as J n (x — > 0) » 
product 7o(P)-/±i (P) contributes to the signal, and we thus need only consider terms with n = 0,m 

with n — 0,m = are DC and do not contribute to the RF spectrum). Also note that, for p <C 1, ^o(P) ~ 1. and /±i(P) 
As an example, for the breathing mode of a silicon nanobeam, one can have C0 o /£2 w 10 s and, with a given by the thermal 
amplitude of oscillation a = ^Jk B T /m efl Q. 2 , a/L M «2x 1CT 13 m/5 x 1(T 6 m= 5 x 1(T 8 , p w 5 x 1(T 3 . Note that this 
condition is independent of the degree of sideband resolution, which is £1/ (k/2); the system can have small modulation in either 
the sideband-resolved or sideband-unresolved regimes. 

In this small modulation approximation, the power oscillating at Q. is 



^ (f )"• In this case, only the 
= ±1 and n = ±l,m = (terms 

=C 2' 
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\ s out,£l\ 



K e Re 



-iQt 



«(P/2> 



-iQt 



-i(£l-A) | + i(-£2-A) 



(E13) 



-/(P/2) e 



(f-fc)(f-i(0- 


A)) 


-/(P/2)e+' a ' 




(|-»A)(| + i(fl4 




/(P/2)e+' a ' 




(f+i(fl-A))(|4 


-zA) 


j(P/2)e-' a ' 





;|-i(n+A))(| + iA) 



This expression can be simplified by combining the terms proportional to (K e /2) 2 (those four terms are really just two terms 
plus their complex conjugates, which can be written as twice the real part of one term). With this simplification, and collecting 
terms proportional to sin(£2f) and cos(£2f ), one finds: 



-= — = cos(£2f) 



^P 

2 



2A 



n-A 



D. + A 



-P 



,2) 

2k 



-A 2 



;|) 2 + (ii-A) 2 (f) 2 + (i2 + A) 2 



(E14) 



-sin(Qf) 



^p 

2 



4 \K 2 +4A 2 J V(f) 2 +(ii-A) 2 (f) 2 + (& + A) : 
k/2 k/2 \ 



+4 



AO. 



,2 f*\ 2 

Q. — A 



(a-A) 2 (f) z + (i2 + A) 2 / 



Cl + A 



4 Vk 2 +4a 2 ; V(f) 2 +(n-A) 2 (f) 2 +(n+A) : 



If we say then, that 



\Sout,a\ 



■■ cos(Q*)A cas + sin(Qt)A sin , 



(E15) 



then the total electrical RF power detected at frequency £2 by a spectrum analyzer (i.e., in both quadratures) is proportional to 

(|wa| 2 ) 2 -(k| 2 ) 2 (AL+A 2 „). (E16) 



4. Optical Forces 



Referring to equations (D3 i and (E8 1, the optical force is equal to 

|2 _ /("-'")y„(p)y m (p) 



2 

\c p (t)\ /L OM = 



i(n—m)Q& 



2Lqm 



n,m \ 2 



-i(nn-A)) (|-i(m£l-A)) 



(E17) 



We will again work in the limit of small modulation depth (p -C 1). In this caase, the DC force is found by considering the term 
for which m = n = (the next DC term is m = n = ±1, which is of order P 2 ), which is simply proportional to the unperturbed 
cavity energy, 
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1 2(k,/k 2 ) 

•foptical.DC = - T T^>Po 

Lom 1 _|_4 



(E18) 



In the limit of small modulation, just as in the consideration of the RF transmission spectrum, there will be four terms at 
frequency i2, 

ct oc 

fWicai.n = F Q cos(Qj) + F I sm(Ql)=FQ— hf> — 



,|2 *«P 

2Lom 
■sin(Q?) 



cos (Of) 



2k 



Q 



Q 



k 2 + 4AV V(i) 2 + (£}-A) 2 (|) 2 + (i2 + A)2 
4£l \ ( Q.-A Q. + A 



K 2 +4A 2 J V(f) 2 + (£i-A) 2 (|) 2 + (n + A) 2 
where Fq and Fj are the in-quadrature and in-phase components of the force, respectively. 



(E19) 



(E20) 



We can now rewrite the differential equation for the acoustic mode, Eq. ( D6 1, to include the effects of the optical mode 



? 1 
a + r,d + iroc= — 



a a 

v Q-pz \-F[ \-Fdc 



Thus, we can rewrite the acoustic mode amplitude's differential equation as 



with 



a+(r i + r)d+(i2 2 + 8i2 2 )a= F DC , 



(E21) 



(E22) 



1 



8(a 2 ) « 2D.8D. 



1 

Fi 

'MeffOCo 



(E23) 
(E24) 



C0 n 



2k, 1st 



k/2 



k/2 



£2L 2 M m eff I K 2 +4A 2 J W^ + ^-A) 2 (f) 2 + (a + A) 2 



8i2 



C0 o 



2kJs| 



a-A 



D. + A 



2i2L 2 M m eff U 2 + 4A 2 y V(f) 2 + (H-A) 2 (f) 2 + (H + A) 2 



(E25) 
(E26) 



5. Power transfer and effective temperature 

The power transfer between the optical and mechanical mode is (see, for instance, |31T ) 



(P) = (/Vicar a) = (ao^m (F e cos(nf) 2 +/7sin(af)cos(af))) 



a £2 n 



■g- 



With equation 



E23 



(E27) 



we can thus write (P) = —^■m e ffD. 2 n T = — (a 2 )m e ffi2 2 r. The free evolution of the average mechanical 
energy, E m obeys 4 (£„,) = —Ti{E m ) + ksTRTi, where Tr is the reservoir temperature. This gives the expected steady-state 
result that (E m ) = I<bTr. In the presence of the optical field, the loss rate of the mechanical system is modified. However, the 
optical field does not modify the reservoir temperature OTI . which means that the evolution of the mechanical energy changes 
to gj (E m ) = — (r, +r) (E m ) +kBTnTi. Thus the steady-state result gives an effective temperature T e ff = jJ^pTs. The effective 
temperature changes the total power of a mechanical mode, whereas the change in linewidth from T just narrows (or broadens) 
the mechanical resonance. 
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6. Calculating the Power Spectral Density 

The thermal amplitude of oscillation is defined by 



\m en a 2 al ermal = l -k B T . (E28) 



Thus we can define a thermal modulation index given by 



P thermal ^0o /r^^n\ 
thermal = ~ T 77 ■ (E29) 



Using equations (E14i - (E16 1 and fithermah one can find the total transmitted power oscillating at frequency v = £l/(2n), \s outy \ 2 . 



Only a fraction of this output power reaches the detector, and we'll call this P@dety = A* |*o«r,v| 2 

The transimpedance gain of the detector, Gti.o, converts optical power to a voltage; the gain has a simple pole, however, at 
frequency Vdet- Thus the total gain of the detector at the frequency of the oscillator, v, is Gti{v) — Gti.q/ (l + (v/v^;) 2 ). This 
voltage is then fed to the input of a buffered channel amplifier with unity gain and a simple pole at frequency, V scope (however, 
in post-processing, the oscilloscope flattens its own response; so that power spectral densities do not contain the oscilloscope's 
pole); then the voltage across the load resistor (input impedance Z) at the output of the channel amplifier is used to compute a 
"power", such that 



p ^=( "lT:f GT i^ ) 2 / z - < e3 °) 

\ (l + (v/v cfa ) 2 ) J 

Because the spectrum of a mechanical resonator is distributed over all frequencies, we must calculate how much power is in a 
given frequency interval, dv = dO/(2jt). If the spectrum is a Lorentzian, we must have that 

P RRV = f s 2 (v)dv' = f ^ (v) x2 dv' = 7tf 2 (v)Sv/2 = . (E31) 

J-°° J -00, fv'-v\ 2Qm 

l+ {5v/2j 

Thus the power spectral density at the peak is 

s2 = 2Q m P RFy 
JTV 

Note that an oscilloscope typically displays S 2 (v) (RBW), where RBW is the resolution bandwidth (the reciprocal of the time 
record length of the FFT). 



7. Extracting the product m e ffL,Q M from experimental RF spectra 

The spectra of FigJ3b of the main text are fit using a sum of Lorentzian lineshapes plus a quadratic background. From the 
fit parameters, Eq. |E32[ allows simple extraction of the total RF power in each mechanical mode. Equations ( |E14[ > - ( |E16| l can 
then be used to calculate the product m e ffLQ M for each mode, albeit indirectly (since it is not simple to invert the equations). 
Given K e , K, A, v m , and v , which are all experimentally-measured parameters (and A is directly set to be the A which gives the 
maximum transduction) the product m e ffLo M uniquely determines the transduced power due to Brownian motion in each mode. 
Thus, by varying m e ffLQ M , one can find the value for the product that gives the correct transduced power. 
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